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ABSTRACT 

We compare spherically-averaged radial number counts of bright (> 10 5 i Q ) dwarf 
satellite galaxies within 400 kpc of the Milky Way (MW) and M31 and find that the 
MW satellites are much more centrally concentrated. Remarkably, the two satellite 
systems are almost identical within the central 100 kpc, while M31 satellites outnumber 
MW satellites by about a factor of four at deprojected distances spanning 100-400 kpc. 
We compare the observed distributions to those predicted for ACDM subhalos using a 
suite of 44 high-resolution ~ 10 12 M Q halo zoom simulations, 22 of which are in pairs 
like the MW and M31. We find that the radial distribution of satellites around M31 
is fairly typical of those predicted for subhalos, while the Milky Way's distribution 
is more centrally concentrated than any of our simulated ACDM halos. One possible 
explanation is that our census of bright (> 10 5 L©) MW dwarf galaxies is significantly 
incomplete beyond ~ 100 kpc of the Sun. If there were ^8 — 20 more bright dwarfs 
orbiting undetected at 100-400 kpc distance, then the Milky Way's radial distribution 
would fall within the range expected from subhalo distributions and also look very 
much like the known M31 system. We use our simulations to demonstrate that there 
is enough area left unexplored by the Sloan Digital Sky Survey and its extensions that 
the discovery of ~ 10 new bright dwarfs is not implausible given the expected range 
of angular anisotropy of subhalos in the sky. 

Key words: Galaxy: halo — galaxies: individual: M31 — galaxies: dwarf — Local 
Group 



1 INTRODUCTION 

The satellite galaxy populations around the Milky Way and 
M31 provide important laboratories for galaxy formation on 
small scales and benchmarks for cosmological p redictions 
assoc i ated with da rk matter halo substructure l|Kravtsovl 
l20ld : lBullockll2010l) . Perhaps the most famous example is 
the Missing Satellites Problem, which points to a mismatch 
in the observed counts of dwarf satellites compared to the 
predicted number of ACDM dark matter sub halos that are 
arguably massive enoug h to have formed stars (jKlypin et al.1 
ll999l ; lMoore et ai1ll999l) . Over the last decade, discoveries of 
new dwarf satellite galaxie s around the Milk y Way in Sloan 
Digital Sky Survey (SDSS: lYork et alJ | 200Ch data haye alle- 
viated this problem to some extent dWillman et al.l 120051: 



Zucker et al. 


20063; IZucker et al 


2006b; 


Belokurov ct al. 


20061: Grillmair 2006; Walsh et alj|2007t llrwin et al.l 12007). 



Specifically, a new population of ultra-faint dwarfs (L ~ 
10 2 — 10 4 Lq) have been discovered that that pile up along 
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the radial detectio n limit of the SDSS (|Kopqsov et alj|200d 
IWalsh et al.l 120091 ) . Completeness corrections suggest that 
there could be as many as a few hundred ult ra-faint dwarfs 
within the virial radius of the Milky Way (|Tollerud et al.l 
120081 : IWalsh et al.ll2009h . 

Given that our Galactic satellite census suffers from 
radially-biased incompleteness effects, M31 and its satellite 
system provides a particularly useful comparison set. The 
observational biases with M31 satellite searches are quali- 
tatively different from those faced in Milky Way. In partic- 
ular, we view M31 from the outside, yet are close enough 
to detect fairly low luminosity dwarf spheroidal galaxies 
there. The Pan-Andromeda A rcheological Survey (PAn- 
dAS; iMcConnachie et al.l 120091) provides a uniform search 
for dwarfs around M31, and has thus far discovered about 
half of the known M31 satellites. Their coverage is com- 
plete to within a projected di stance of 150 kpc and outer 
M31 halo coverage is ongoing (jRichardson et ahl 1201 ll ). In- 
terestingly, the current M31 dwarf count, which is almost 
certainly incomplete beyond 150 kpc, is larger than the 
Milky Way count by about a factor of two at fixed lumi- 
nosity (see below). This difference is striking given that 
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the two spirals have the sa me luminosity to within ~ 30% 
(|van den Bergh et alj|200ch . 

Beyond the question of overall counts, another useful 
comparison concerns the relative radial distribution of satel- 
lite galaxies. The Milky Way satellites are known to be more 
centrally concentrated around the Galaxy than would be 
expected for subhalos in dissipationless ACDM simul ations 
(|Moordl200il ; IWillman et alJl2004lMacci6 et alj|20ich . This 
basic mismatch has been used as a potential testing ground 
for ideas to solve the Missing Satellites Problem. For exam- 
ple, early-forming halos that formed before reionization tend 
to be more centrally con c entrated than the overall popula- 
tion (|Bullock et al.ll2000l ; iMoord 120011 ). The problem with 
this explanation is that the mismatch between naive theory 
and observation ap p ears t o be much less severe for M31. 

IWillman et~atl l|2004l were the first to emphasize that 
known M31 satellite galaxy population was somewhat more 
extended than that of the Milky Way, and suggested that the 
disagreement could be explained by incompleteness in the 
outer Milky Way halo, possibly as a result of the difficulty 
in detecting low-surface brightness systems at low Galactic 
latitude. Importantly, at the time of their work, all of the 
known satellites of M31 and the Milky Way were "bright" by 
today's standards, with luminosities L > 10 5 L©. They em- 
phasized that this possible incompleteness in (bright) satel- 
lites provided motivation for searches for new dwarfs via 
resolved-stars in the Sloan Digital Sky Survey. As already 
mentioned, these searches have proven immensely success- 
ful at finding new faint dwarfs, but new bright dwarfs of 
the type known before 2004 have remained sparse: all but 
one of the ~ 15 new dwarfs that have been discovered since 
2004 are less luminous than 10 s Lq and none are more lu- 
minous than 3 x 10 Lq. Over the same period, 20 new M31 
satellite galaxies have been discovered, 14 of which are more 
luminous than 10 Lq . 

In this paper we reexamine the radial distributions of 
M31 satellites and Milky Way satellites. We restrict our- 
selves to the brightest systems: bright enough that we could 
possibly assume our census of them complete within the 
Milky Way virial radius. Given the new discoveries of mainly 
bright satellites around M31, we expect to see the dif- 
ference in their radial populations t o be even more dis - 
crepant than previously quantified (|McConnachid 120061 ). 
This discrepancy may grow yet with the discovery of ad- 
ditional satellites aroun d M31 as PAndAS coverage expands 
(jRichardson et alj|20lj ). 

The outline of this paper is as follows: in Section 2 we 
describe the data sets of the MW, M31 and our suite of 
ACDM simulations. Section 3 contains the results of our 
comparison, Section 4 discusses the possible effects of in- 
completeness in the satellite populations and in Section 5 
we conclude. 



2 DATA 

2.1 Satellites of the MW and M31 

The M31 and Milky Way satellite data sets we examine in 
this paper are summarized in Table 1, with data taken from 
the review bv iMcConnadfil Galaxies are listed in 



the order they were found in order to emphasize the great 
explosion of discoveries in the post SDSS /PAndAS era (hori- 
zontal break). We include only systems that could be within 
400 kpc physical separation of each host, given published 
distance uncertainties. 

The Milky Way satellites listed in Table 1 are restricted 
to those brighter than the faintest known M31 dwarfs (> 
10 4 Lq), 15 in all. Twelve of these are brighter than 10 5 Lq, 
the approximate completeness limit of SD SS at 400 kpc and 
of PA ndAS for M31 satellite discovery l|Richardson et all 
l201lh ■ Only one of the new dwarfs discovered in the ~ 1/3 
of the sky covered by the SDSS is brighter than this: Canes 
Venatici I at 220 kpc distance. This has provided some mo- 
tivation for assuming that the Milky Way census of bright 
dwarfs is reasonably complete out to the virial extent of 
the Milky Way (~ 300 kpc), though it will be impossible 
to quantify this completeness with certainty until uniform, 
full-sky, resolved-star searches have been completed. 

The M31 data set includes a total of 32 galaxies, more 
than half of which were discovered in the last six years. The 
census of dwarfs in the vicinity of M31 is well-understood 
withi n a projected radius of about 150 kpc of the M31 cen- 
ter l|Richardson et all 1201 ll ). Although the faintest known 
M31 dwarfs known have luminosities down to ~ 10 4 Lq, 
the PAndAS surv ey's limit for efficient d etection is closer 
to L v ~ 1O 5 L l|Richardson et al J [20111 ). This complete- 
ness limit applies only to the region within 150 kpc of the 
M31 disk. Of the 32 M31 satellites listed, 26 are brighter 
than 10 5 Lq. We will use this luminosity as a characteris- 
tic luminosity for comparison in this work. We have not 
included And VIII a nd And XVIII in our comparison set. 
iMerrett et al.l (|2006l ) have questioned And VIII's classifica- 
tion as a galaxy because of tidal features. And XVIII has 
a very high deprojected separation from the center of M31, 
falling outside of our 400 kpc region of consideration. 

The deprojected distances and errors in Table 1 were 
derived from the hel iocentric distances and errors quoted in 
iMcConnachid |2012i ). The separation between M31 and one 
of its satellites is given by 

tmzi = VR 2 + d 2 - 2 Rdcos 6 , (1) 

where R is the heliocentric distance of M31, d is the heliocen- 
tric distance of the satellite and 6 is the angular separation 
between the two. We sampled the heliocentric distance of 
M31 and of each satellite 1000 times, drawing from a Gaus- 
sian based on quoted line-of-sight errors to obtain the error 
limits of r_M3i, which are not always Gaussian-like them- 
selves. 

2.2 Simulations 

We use a total of 44 Milky-Way size halos simulated as 
part of the ELVIS project, which is described in detail in 
Garrison-Kimmel et al. (in preparation); here we summa- 
rize the relevant properties. ELVIS is a suite of collisionless 
zoom-in simulations designed to study the Local Group. It 
consists of twenty-two halos in paired systems that were cho- 
sen to resemble the M31 and Milky Way in mass and phase- 
space configuration in addition to twenty-two halos that are 
isolated, mass-matched analogues, which serve as a control 
sample to examine any trends due to the paired nature of 
the Local Group system. The host dark matter halos have 
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Table 1. Satellite galaxies of M31 (left) and Milky Way (right) used in this work, listed in order of their discovery date. The horizontal 
band marks the transition from the "classical" era of discovery to the post-SDSS/PAndAS era. With a deprojected distance greater than 
400 kpc (with relatively small uncertainty), And XVIII may not be a bound satellite of M31; it is therefore not listed here or included in 
any figures. For completeness, we list all satellites brighter than Lv = lO 4 Z/0 here, though we only present data for systems fainter than 
Ly = 1O 5 L0 i n Figure [T] We use th e mor e conservative cut at Ly = lO 5 Z/0 for subsequent figures and radial comparisons. Luminosities 
are taken from lWatkins et al. | J2012h and lMcConnachiei J2012T) while the host-satellite separations and associated errors were computed 
from the heliocentric distances quoted by those authors. 



virial masses between 1.0 — 2.8 x 1O 12 M0 and associated 
virial radii i? v i r = 265 — 370 kpc, with maximum circular 
velocities that range from Knax = 155 — 225 kms -1 . 

All ELVIS halos were simulated in a ACD M cosmology 
with parameters based on WMAP7 results (|Larson et al.l 
1201 lb : a s = 0.801, fl m = 0.266, ft A = 0.734, n 3 = 0.963, 
and h = 0.71. The simulations were run with identical mass 
resolution (particle mass m p — 1.9 x 10 5 Mq); likewise, all 
simulations were run with a Plummer-equivalent force soft- 
ening of 140 pc, which was comoving until z = 9, after which 
a physical softening length was imposed. All simulations are 
uncontaminated by low-resolution particles to at least 900 



kpc from the center of e ach host. Halo substruc ture was iden- 
tified with Rockstar IjBehroozi et al.l l2012a[). and merger 
trees w ere built using consistent-trees (jBehroozi et al.l 
l2012bh . At this resolution, we are complete to subhalos with 
current maximum circular velocity V max > Skms^ 1 and to 
those with Vp Ca k > 10 kms" 1 , where VJ, e ak is the largest cir- 
cular velocity ever achieved by the main branch of a subhalo 
progenitor. 

Additionally, two of our host isolated halos were resim- 
ulated at an even higher resolution with a particle mass of 
m p — 2.35 x 10 4 Afo and a force resolution of 70 pc, com- 
parable to the resolution of the level 2 Aquarius simula- 
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Figure 1. Satellite galaxy luminosity functions of the MW (blue) and M31 (red) for all galaxies within 400 kpc of either host (left); 
those within 100 kpc of either host (middle); and those in the radial range 100-400 kpc (right). While the luminosity functions within 
100 kpc are very similar, M31 has approximately four times as many satellites at fixed luminosity in the 100-400 kpc range. 



tions (Springcl ct al. 2008) and the Via Lactea I simulation 
(|Diemand et al.l 20071) . The velocity functions and radial dis- 
tributions of the two higher resolution host halos are nearly 
identical to those of the lower resolution counterparts to 
Vpoak values well below those used for selecting subhalos in 
this work. 

In our fiducial radial distribution comparisons, we se- 
lect the 30 subhalos with the highest V^ ea k within 400 kpc 
of each host. This choice is motivated by the fact that M31 
has approximately 30 satellite galaxies within this radius; 
this is also the approximate distance to which SDSS is com- 
plete to L > 10 5 Lq dwarfs. In practice, this amounts to 
selecting halos with V^ ea k values greater than 20kms _1 to 
33kms _1 , depending on the host. The overall shapes of the 
distributions are fairly insensitive to this specific choice. For 
example, when we chose the 100 most massive subhalos, or 
when we use a fixed Vp Ca k cut, we find that the radial dis- 
tributions of subhalos do not change significantly. In the 
Appendix, we explore the spatial distributions of subhalos 
with the 30 largest values of Knax and those with the 30 
largest values of z pca k (the redshift at which Vpeak is at- 
tained) . For the latter case, we restrict ourselves to subhalos 
with V m ax > lOkrns -1 . 



3 RESULTS 

3.1 The MW and M31 satellite systems compared 

The satellite luminosity functions of the MW and M31 are 
shown in Figure [1] The left panel includes all satellites 
within a deprojected radius of 400 kpc of either hostQ Con- 
sidering this full sample, M31 has approximately a factor of 
two more satellites at fixed luminosity. Intriguingly, however, 
the two systems have luminosity functions that are essen- 
tially identical when only satellites within 100 kpc are con- 
sidered (middle panel). The difference in satellite luminosity 
functions between the two systems is entirely attributable 
to differences at large radial separations (100 - 400 kpc), as 



1 In this figure we have used the best-fit radial distances of satel- 
lites listed in Table 1. 



shown in the right panel of Figure [T] While M31 has ap- 
proximately 20 satellites brighter than 10 5 Lq in this radial 
range, the Milky Way has only 5. We note, however, that 
the shape of the luminosity function remains very similar 
at these large radii, even though the amplitude differs by a 
factor of ~ 4. 

A complementary demonstration of the difference be- 
tween the satellite systems of the two Local Group giants 
is given in Figure [2] which compares the cumulative radial 
count of satellites brighter than 10 Lq around each host. 
Solid lines show median distributions and the dashed lines 
show 68% uncertainties derived using Monte Carlo realiza- 
tions of the deprojected distances and errors listed in Table 
1. M31's satellite population is significantly more extended 
than that of the MW, almost entirely due to differences at 
> 100 kpc. 

In the context of observational completeness, it is 
remarkable that M31 has so many more satellites at a 
large galactocentric distances. PAndAS has provided a 
uniform search for (bright) dwarf galaxies only within a 
projected distance of 1 50 kpc of the Andromeda galaxy 
(jRichardson et al.ll201ll ). so it would not be surprising if 
there are even more M31 dwarfs waiting to be discovered 
beyond 150 kpc - a region marked by the vertical dashed 
line in Figure [2] Any such discoveries would only enhance 
the discrepancy between the MW and M31. 

One way the difference in cumulative radial distribu- 
tions might be reconciled is if our M31 census were incom- 
plete at small projected radius, perhaps because of the diffi- 
culty in detecting low surface brightness dwarfs in the vicin- 
ity of M31's disk. Interestingly, of the M31 dwarfs less lu- 
minous than 10 6 Lq, And IX is the closest at a projected 
distance of R = 36 kpc. Meanwhile, there are 2 (or 3, if And 
VIII is included) very luminous satellites that are within this 
radius. If there were a few more low-luminosity dwarfs hid- 
ing in glare of the M31 disk, this would act to shift the M31 
satellite distribution towards a more concentrated shape, 
more in line with that of the MW. However, discoveries 
of this kind would further exacerbate the overall difference 
in normalization between the two satellite luminosity func- 
tions. Moreover, additional small-R discoveries around M31 
would also drive the M31 radial distribution away from the 
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Figure 2. Cumulative number of satellites with L > 1O 5 L0 
within a deprojected distance r from either the MW or M31. The 
solid and dashed lines show the median and 68% confidence inter- 
vals distributions implied by current observational uncertainties. 
The spatial distributions are similar within ~ 100 kpc, while the 
M31 satellite system is significantly more extended than that of 
the MW outside of this radius. This is in spite of the fact that 
the M31 system is almost certainly incomplete at distances on 
the right of the vertical dashed line, which marks the edge of the 
region within which the PAndAS project has provided a uniform 
search for M31 dwarfs. 



expectations for subhalos derived from our ACDM simula- 
tions. As we demonstrate in the next section, the distribu- 
tion of presently-known M31 satellites is much more in line 
with subhalo distributions derived from our LCDM simula- 
tions than the known satellite distribution around the Milky 
Way. 



3.2 Comparison to Simulations 

Since the Milky Way and M31 have such different radial dis- 
tributions of satellites, it is natural to ask which (if either) 
is more in line with theoretical expectations. We expect a 
correspondence between sub halos with large values of Vpea k 
and luminous satellites (e.g., Bullock 2010; Kravtsovll2010l ). 
Figure [3] therefore compares the observed radial distribu- 
tions of satellites of the MW (blue) and M31 (red) to those 
of the 30 top Vp Ca k subhalos from our 44 ACDM simulation 
(gray lines); for each data set, the radial profile is normal- 
ized to the cumulative count within 100 kpc. As in Figure(2] 
the dashed colored lines show uncertainties associated with 
deprojection for the observational data. 

Although there is significant scatter in the number of 
satellites at large radii (> 300 kpc) in the simulations, it 
is clear that the M31 system is consistent with the profiles 
expected from subhalo distributions while the MW is more 
centrally concentrated than any of the simulated systems. 
A typical simulated halo has 4-5 times as many satellites 



Figure 3. The cumulative number of satellites within radius r, 
normalized to the count within r = 100 kpc separation. The solid 
(dashed) blue and red lines show the median (68% spread) profiles 
for the MW and M31 satellite systems while the grey and black 
lines show the profiles of subhalo distributions derived from 44 
ACDM simulations. The subhalo distributions were made from 
the 30 subhalos with the highest V pca ]j from each simulation. 
The Milky Way's distribution is more centrally concentrated than 
any of the simulated systems, while M31 looks fairly typical. For 
reference, the five solid black lines show subhalo profiles for the 
five highest virial mass halos in our sample (2.7 - 2.8 X 10 12 M Q ) 
and the black dotted lines lines are from our five lowest mass 
halos (1.0 - 1.2 X 10 12 Af Q ). While there is some tendency for the 
highest mass systems to have more extended distributions (when 
measured this way) , the correlation is weak and shows substantial 
scatter. 

within 300 kpc as within 100 kpc, while the corresponding 
number for the Milky Way is only 1.5. 

The solid (dotted) black lines in the Figure correspond 
to the five highest (lowest) mass host halos in our sample. 
While two of the five most massive halos also have subhalo 
populations that are among the most extended of our sam- 
ple, as might be expected based on their large virial radii, 
three of the five most massive hosts have subhalo popula- 
tions that are either fairly typical or even somewhat more 
concentrated than average. Any trend with host halo mass 
is fairly weak and is overwhelmed by sc atter at fixed mass 
(this was also seen bv lWang et ai]|2013| y 

It is useful to have a parameterization of the degree of 
concentration of a halo's satellite population. For this we 
define a parameter B400 to be the ratio of the number of 
satellites within 400 kpc to the number within 100 kpc: 

B400 = N(< 400kpc)/Af(< 100 kpc). (2) 

B400 is then a measure of the "puffmess" of the subhalo 
population. The distribution of B400 values for all of our 
simulated host systems is shown by the solid histogram in 
Figure [J] while the measured values for the MW and M31 
are indicated by arrows. As expected, we see that the satel- 
lite distribution of the MW is more concentrated than any 
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Figure 4. The distribution of subhalo radial distribution shapes 
- B 40 o = N(< 400kpc)/Af(< lOOkpc) — for all of our simu- 
lated halos (solid line) with the observed values for the MW and 
M31 indicated by arrows. The shaded histogram shows the subset 
of simulated halos that are in Local Group like pairs. The MW 
is clearly an outlier from the simulated distributions in that its 
satellite galaxy population is much more concentrated; M31 is 
more typical of the simulated systems. 



of the simulated subhalo systems, while M31 is fairly typ- 
ical of what is seen in the simulations. The shaded blue 
histogram shows the subset of our simulated halos that are 
paired, Local Group analogs; these halos are fairly similar 
to the isolated, mass-matched analogs in their B400 distribu- 
tions. As a result, the fact that M31 and the MW are paired 
should not significantly skew their satellite radial distribu- 
tions compared to field systems (at least when spherically 
averaged; we will explore angular anisotropy in more detail 
in a future paper). 

Even though the MW is an outlier relative to all of our 
simulated systems, one might posit that the high concentra- 
tion of the Milky Way's satellite population is attributable to 
some physical process. However, we were unable to find any 
correlation between -B400 and halo concentration, mass ac- 
cretion history, or virial mass. In the Appendix, we show ex- 
plicitly that the radial subhalo distributions are not tightly 
correlated with the bulk infall time distributions of satellites. 
Given that we are unable to provide a natural explanation 
for why the MW system is highly anomalous compared to 
predicted subhalo distributions (while M31 is not), in the 
next section we explore the possibility that the census of 
bright MW satellites may be incomplete at radii beyond 
100 kpc. 

In the Appendix we explore a few other options for 
selecting subhalos and the associated implications for ra- 
dial distribution. For most reasonable choices, the selection 
on Vpeak we use for our fiducial analysis produces a more 
centrally concentrated subhalo distribution than those with 
the highest Vmax- However, we show that by sub-selecting 
Knax > lOkms -1 subhalos that reach their Vp ea k values the 



earliest (i.e. were accreted into the host halo the earliest) 
we can produce subhalo di stributions that a re more cen- 
trally concentrated (see also lWang et al.ll2013T ). and more in 
line with the observed radial distribution of the Milky Way. 
However, this model is not very well motivated as a physical 
selector for luminous subhalos. As shown in the Appendix, 
it effectively picks out ~ 10% of the subhalos that have 
Vpcak < 30km s" 1 , and ignores an equal number of systems 
with Vpoak > 30km s _1 . There is no astrophysical reason to 
expect these most massive progenitors to be dark. 

3.3 Are there undiscovered bright satellites 
around the Milky Way? 

The Milky Way satellite census may still be incomplete at 
large radii, even to satellites similar in luminosity to the 
classical dwarfs (L > 10 5 Lq). The only truly uniform survey 
with well- understood completeness limits is SDSS/SEGUE, 
which has covered approximately 1/3 of the sky. One bright 
dwarf (L > 10 5 Lp), Canes Venatici I, was discovered in 
this area (|Zucker et alj|2006b[) . Is it therefore worthwhile to 
ask the question: how many new Milky Way dwarfs could 
plausibly be undiscovered? We investigate this question first 
in terms of sky coverage and angular anisotropy, then in 
terms of radial distributions. 

Of the bright MW satellites, there are 6 that lie in- 
side of the SDSS footprint. We can be confident that our 
census is complete within this region. Given this, we can 
provide a rough estimate for a plausible number of bright 
dwarfs that could exist outside this footprint using our sim- 
ulations. Specifically, for each of our 44 systems, we count 
subhalos in randomly-oriented angular cones the size of the 
DR8 footprint (14555 deg 2 ) with a depth of 400 kpc. We 
determine a Vp ea k cut that yields 6 subhalos within this re- 
gion. We then determine the total number of subhalos over 
the whole sky with Vpeak higher than the cut. This gives the 
total expected number of bright satellites within a spherical 
region of 400 kpc, iVexpected, based only on the count within 
an SDSS area. 

For each host halo, we perform 10 4 random realizations 
of mock pointings to derive an A^xpected distribution based 
on our full sample. The results are shown in Figure The 
median prediction is 18 bright satellites over the whole sky. 
This can be compared to the 12 known bright MW satellites, 
so it is not implausible to think that ~ 6 undiscovered bright 
satellites of the Milky Way could be out there. In fact, the 
distribution extends to at total ~ 40 satellites, which would 
correspond to a total of ~ 28 undiscovered satellites. 

How many satellites would need to be discovered be- 
yond 100 kpc if the MW satellite distribution were to match 
the expectations of our ACDM simulations? The Milky Way 
has 7 bright satellites within 100 kpc, so for each of our sim- 
ulated halos, the value of -B400 can be used to compute the 
total number of bright MW satellites expected based on the 
known number within 100 kpc: 

A'expcctod = N M w(< 100 kpc) x B400 = 7B 4 oo • (3) 

The right panel of Figure [5] shows that this predicts that 
there would need to be a total of 20 to 40 bright MW satel- 
lites, substantially more than the 12 that have been discov- 
ered so far. The implication is that there would need to be 
between ~ 8 — 28 new satellites beyond 100 kpc in order for 
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Figure 5. Left: The total number of Milky Way satellites expected over the whole sky given the observation that there are 6 bright 
satellites in the SDSS footprint (14555 deg 2 ). The calculation uses subhalos in our ACDM simulations to correct for areal incompleteness, 
allowing for anisotropy. The assumption is that only the SDSS area is known with confidence to be compete to systems with L ~ 10 5 Lq. 
Right: The total number of satellite galaxies expected with 400 kpc of the M W, estimated from the number of MW satellites known 
within 100 kpc and then assuming that ratio of satellites within 400 kpc to those with 100 kpc matches the the expectation for ACDM 
subhalos (see text). The shaded regions represent the contribution due to the paired "Local Group-like" simulations. The blue arrows 
indicate the number of known MW satellites brighter than 10 5 Lq. 



the Milky Way to fall within the range of our simulations. 
Interestingly, this number overlaps with the allowed range 
from the angular coverage correction discussed above (and 
in the left panel of Figure [5]) . 

To have escaped previous detection, a bright satel- 
lite would likely need to have a large half-light radius 
(~ 1 kpc) and low surface brightness. Canes Venatici I, 
the bright satellite discovered in SDSS, does indeed have 
such characteristics: its central surface brightness is /i = 
27. 1 mag/arcsec 2 and its projected half ligh t radius is 560 
pc jZucker et al.ll2006bl ; lMcConnachiell2012r i. 



4 DISCUSSION AND CONCLUSIONS 

The results of Section[3]show that the satellites of Milky Way 
and M31 have very different radial distributions. Within 100 
kpc of their respective hosts, the L > 10 5 L© satellite galaxy 
populations of the two galaxies have similar radial distribu- 
tions and luminosity functions, whereas between 100 and 400 
kpc, M31 has approximately four times as many satellites as 
the Milky Way. While ACDM simulations predict substan- 
tial scatter in the radial distribution of subhalos, the MW's 
satellite population is more concentrated than all of the sim- 
ulations we have investigated. This result holds for reason- 
able associations between subhalos and galaxies, though ex- 
treme subsets of subhalos, chosen on formation time, can 
alter this conclusion (see below). M31, on the other hand, 
looks fairly typical when compared to the simulations. 

It is possible that the Milky Way is simply an outlier in 
terms of the spatial distribution of its satellites - though it 
would need to be extreme enough to be rarer than one out 
of the 44 systems we have simulated. A more intriguing pos- 
sibility is that the census of Milky Way satellites is incom- 
plete not just for the faintest dwarfs, but also for brighter 
systems (L ~ 10 5 L@). The discovery of Canes Venatici I by 



IZucker et ail (|2006bh in SDSS data indicates that this is in- 
deed a distinct possibility. Using the distribution of subhalos 
from our simulations, we might expect 8-20 additional satel- 
lites if the radial distribution of the MW satellites is similar 
to those of subhalos in our simulations. We showed that the 
angular anisotropy of subhalos in our simulations allow for 
the possibility of 8-20 additional satellites in the area of the 
sky not yet surveyed for new dwarfs by the SDSS. 

If the shape of the radial distributions of subhalos is 
a template for the radial distribution of satellite galaxies 
about their host, then the Milky Way almost certainly has 
new bright satellites yet to be discovered. This would be 
consiste nt with the expe c tation s of hydrodynamical simula- 
tions of lBovill fc Ricottl (|201ll ) who have argued that there 
should be more bright satellites around the Milky Way based 
on their models. While the simulations we have used for our 
comparison are dissipationless, including baryonic physics is 
likely to exacerbate the discrepancy between the predicted 
and observed distribution of satellites in the Milky Way, 
as the MW disk can d eplete substructure that pass nearby 
(|D'Onghiaet alj|20ld ). 

One way to increase the predicted concentration of the 
MW satellite population, and therefore to bring the MW 
into better agreement with dissipationless ACDM simula- 
tions, is to associate the Milky Way's bright satellites with 
its earliest forming subhalos having Vp C ak > lOkrns -1 . Such 
an association also predicts that the Milky Way's bright 
dwarf spheroidals (1) formed in halos with Vp Ca k < 30kms _1 
and (2) reside in halos that currently have Kn ax < 20 kms -1 . 
This sce nario would be in agreemen t with the values de- 
rived by iBovlan-Kolchin et al.l |201ll ) based on kinematics 
of the MW's bright satellites. In this case, there would be 
20-40 dark matter satellites of the Milky Way that are more 
massive but less luminous than its bright dwarf spheroidal 
galaxies. This scenario is highly unnatural in current models 
of galaxy formation. 
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Although we do not see a significant difference in the 
subhalo radial distributions between paired and non-paired 
host halos, we can not rule out the possibility that the evo- 
lutions of the Milky Way and M31 are coupled in some 
way. The asymme trical 3D distribution of M31's satellites 
()Conn et al.ll2013h may hint at such a connection. We are 
currently using the ELVIS suite to evaluate the frequency of 
such "great plane" configurations in ACDM simulations. 
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Figure 6. Subhalo concentration distributions (-B400) plotted as 
a function of 2 pea k> the median redshift at which V pea k is attained 
for each subhalo system. This is analogous to comparing -B400 to 
the median infall time of each subhalo set. The error bars show 
the la spread around the median value of z pea k f° r each sub- 
halo set. There is no evident correlation, and we are also unable 
to find a correlation with any property of the main halo (e.g., 
concentration, mass, or formation time.) 



Figure 7. A histogram of -B400 values for all 44 host halos com- 
puted by selecting the subhalos with the highest Vp ea k (black), 
Vmax (magenta) and 2 pea k (cyan). The range and distribution 
of B400 depend somewhat on how subhalos are selected. Our se- 
lection of subhalos with the 30 highest V pe ak predicts somewhat 
more concentrated distributions of satellites than does the se- 
lection based on Vmax. The most concentrated distributions are 
found when selecting based on the earliest forming (highest z pea k) 
subhalos. 



APPENDIX 

It is reasonable to expect the concentration factor, B400, to 
be correlated to with some property of host halo. For exam- 
ple, the concentrations of dark matter halos are known to 
be correlated with formation redshifts (N avarro et al.|[l997l ; 
iBullock et al.ll200ll ). We have searched for relationships be- 
tween -B400 and several properties of the host halo, including 
virial mass, virial radius, J? ma x, Vmax, and M(< -Rmax), but 
were unable to find any correlations. Figure [S] illustrates one 
such null result: B400 is plotted as a function of the median 
(asterisks) and 68% spread (error bars) in a pGa k, the expan- 
sion factor at which a subhalo reaches V pea k, for each ELVIS 
halo. While the spread in a pea k is large for each system, there 
is no correlation with -8400- 

Figure [7] explores a few different ideas for associating 
subhalos with bright satellites and the effect of that selection 
on their radial distributions. Specifically, for each of our 44 
hosts we plot the implied B400 (puffiness) parameter (Equa- 
tion 2) derived for three different choices of subhalo selec- 
tion: the thirty highest V pca k (black, fiducial); the 30 highest 
Vmax (magenta); and the 30 highest z pca k (cyan; defined as 
the redshift where V pca k was reached). We see that by select- 
ing on z pca k we produce distributions that are not as puffy as 
our fiducial choice, and marginally consistent with the B400 
of the MW, which is 1.7. However, subhalos chosen in this 
manner are an odd subset in terms of their mass. Figure [8] 
shows the V pea k and V max mass functions of one example sys- 
tem selected in this manner (dashed) compared to all of the 
subhalos associated with this host (solid). It would demand 
a large number of V pca k > 30km s -1 and V max > 20km s -1 



subhalos be dark. Such a model would be very difficult to 
understand from a galaxy formation perspective, as there is 
no obvious mechanism for allowing galaxies to form in ha- 
los of 30kms _1 while suppressing galaxy formation in halos 
of 30 — 50kms -1 . Moreover, the choice amounts to select- 
ing subhalos that reached their peak V max value at redshifts 
between 2 and 4, well later than the epoch of reionization, 
which might offer more obvious channels for galaxy suppres- 
sion mechanisms. 
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Figure 8. Cumulative counts of subhalos as a function of V ma x (left) and Vpcak (right) in a representative simulated halo. The solid lines 
correspond to all subhalos within 400 kpc of host halo, while the dashed lines represent a subset of those halos with the highest values of 
^peaki the redshift at which the subhalo reached the Vpeak- The ~ 10 most massive subhalos have V pca ]j > 40kms _1 or V m ax > 30kms _1 , 
whereas the ~10 largest early forming subhalos have Vp ea k > 18kms — 1 or V max > 12kms . 
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